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^ ABSTRACT 

O ' 

■ We analyze the metallicity distributions of globular clusters belonging to 28 early-type galaxies in the 
survey of Kundu & Whitmore (2001). A Monte Carlo algorithm which simulates the chemical evolution 
of galaxies that grow hierarchically via dissipationless mergers is used to determine the most probable 

^ ' protogalactic mass function for each galaxy. Contrary to the claims of Kundu & Whitmore, we find that 

OO , the observed metallicity distributions are in close agreement with the predictions of such hierarchical 

■ formation models. The mass spectrum of protogalactic fragments for the galaxies in our sample has a 
' power-law behavior, n{Ai) oc with an index of a ~ —2. This spectrum is indistinguishable from 

the mass spectrum of dark matter halos predicted by cold dark matter models for structure formation. 
We argue that these protogalactic fragments — the likely sites of globular cluster formation in the 
early universe — are the disrupted remains of the "missing" satellite galaxies predicted by cold dark 
matter models. Our findings suggest that the solution to the missing satellite problem is through the 
' suppression of gas accretion in low-mass halos after reionization, or via self-interacting dark matter, and 

argue against models with suppressed small-scale power or warm dark matter. 

O ' Subject headings: Cosmology: observations, Cosmology: dark matter, Galaxies: Elliptical and 

^ , Lenticular, cD, Galaxies: Formation, Galaxies: Star Clusters 

a 



1. INTRODUCTION problem have focusscd on the properties of the dark mat- 

<^ ij J 1 4-4. //^7-lHJ^ J 1 4.4. 4.- 4.1, ter itself. Competing models suggest that the dark matter 

Cold dark matter (CDM) models otter an attractive the- . / •' a • » 



. ^ ory for the formation of cosmological structure. These warm (Cohn, Avila-Reese & Valenzuela 2000), decay 
U models, and particularly the ACDM variant, show impres- '""^ (p^^ 2001), self-mteractmg (Dave et al. 2001), repuL 
P. ; sive self-consistency with many of the most fundamental "^^^ (Goodman 2000) or annihilating (Riotto & Tkachev 
astrophysical observations: i.e., the expansion age of the 2000). Ahernatively, the standard CDM picture may be 
universe, the primordial element abundance, the two-point modified to include a broken scale mvariance which sup- 
correlation function of cluster and galaxies, the local space P'^esses smaU-scale power (Kamionkowski & Liddle 2000). 
density and rcdshift evolution of massive clusters, and the ^et another modification has been proposed by Bullock 
large-scale streaming motions of galaxies in the local uni- ^t al. (2000), who note that a photo-iomzmg background 
verse {e.g., Bahcah et al. 1999). Nevertheless, it has be- expected to mhibit the accretion of gas onto low-mass 
come clear that this success on large scales {i.e., R ;> 0.1-1 ^alos after the epoch of reionization. In this latter sce- 

T,f \ , , 4.jiiui 1 c 1- f ^ nario, small-scale power is not suppressed, but /izrfden: as a 

Mpc) does not extend to the level ot galaxies (see, e.g.. . , i ,• r t , i i • • i i 

c 1 1 J p 1 onni r ■ \ a- i I 4- j ' mixed population ot disrupted and survivmg dark-matter- 

bellwood & Kosowsky 2001 tor a review). Simply stated, , • , , , i , 

/^TMv/T J 1 J u J f 11 dominated subhalos. 

CDM models produce an overabundance or power on small t i ; <• n . • , . i • r 

, , T J 4. 4. 1 4.1 u J 1 1 In what lollows, we present a simple technique tor mea- 

scales when normalized to match the observed large-scale . n t . i ii i • i 

J , -1111 J -i surmg the mass distribution 01 disrupted subhalos m low- 

power: e.g., the models yield dark matter density pro- , ,"n , . ^ , i • ,i 

ci 4.1 4. 4. ii, 4-1, u J /T71 on- redshitt galaxies. Our technique relies on the GC systems 
nles that are steeper than those observed t lores & Pri- °, . ,., , r .i ■ . 

1 inn /I \ J J- i 4-u u -c J r 1 • ot thcsc galaxies which, by virtue oi their extreme age 

mack 1994), and overpredict the number 01 dwart galaxies , " ^, . ' , ^ ^. ^ . '^^ 

rr Tin,-4- p /-I -J • 1nno^ inr • 4-1, and apparently universal formation cthciency, represent 

(Kauttmann, White & Guiderom 1993). We examine the ■ , , , r .^ ■ ■ ^ i r i i 

, . 11 J a • • 4- 11-4- ui » ideal tracers ot the original population ot low-mass sub- 

ter issue — the so-called missmg sateUite problem , , . , . , , , ; , ;• r 

, • 4-1, 1 1, 1 14- • 4- J -4-1, halos or, m the terminology used here, protogalactic trag- 

by using the globular clusters (GCs) associated with a , ' , , , ; , • , r i 

ments. ihe observational database consists ot GC color 



, , , . ' 11 J a • • 4- 11-4- ui » ideal tracers of the original population of low-mass sub- 
latter issue — the so-called missmg satellite problem , , . , . ° , , . , r 

, • 4-1, 1 1, 1 14- • 4- J -4-1, halos or, m the terminology used here, protogalactic trag 

— by using the globular clusters (GCs) associated with a , ' , , , ; , • , r i 

, r 1 1 1 • i J • 4-1, ments. ihe observational database consists ot GC coloi 

sample ot nearby, early-type galaxies to derive the mass , n- i- ^ -i r on i ^ i • 

J.,.,,. r 4-i ^-i-4-i i and metallicity distributions lor 28 early-type galaxies pre- 

distribution or accreted material m these systems. ^ 1 Txn /onniN Al/l 1 

,j , , , , 1 • ii • • X 11-4- sented by Kundu & Whitmore (2001). Although most cos- 

Most recent attempts to explain the missmg satellite \ / o 
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mological evidence appears to favor the notion that the 
GC systems of large; galaxies were assembled hierarchi- 
cally {e.g., see Gnedin, Lahav & Rees 2001 and references 
therein), Kundu & Whitmore (2001) claim that the GC 
metallicity distributions of their program galaxies are in- 
compatible with such scenarios. In this paper, we analyze 
the data of Kundu & Whitmore (2001) using a quanti- 
tative numerical technique, described in a pair of earlier 
papers (Cote, Marzke & West 1998; Cote et al. 2001), 
which is designed to recover the protogalactic mass spec- 
trum of a galaxy. The method assumes that the mergers 
and accretions by which cause galaxies to grow hierarchi- 
cally induces little or no gas dissipation. The mass spectra 
which we derive for the protogalactic fragments in each of 
these galaxies are indistinguishable from those predicted 
by CDM models for the population of low-mass dark mat- 
ter halos, a finding which has obvious implications for the 
missing satellite problem. 

2. GLOBULAR CLUSTER DATA 

Since GCs as a class are composed of coeval stars, the in- 
terpretation of their broadband colors is far more straight- 
forward than is the case for galaxies. Moreover, for com- 
posite stellar populations with mean ages of a few Gyr or 
more, broadband colors are primarily sensitive to metal- 
licity, not age (Bruzual & Chariot 1993; Worthey 1994). 
Luckily, the great majority of the GCs in these systems 
appear to be older than this, according to optical spec- 
troscopy, color-magnitude diagrams, and luminosity func- 
tion turnover magnitudes for GCs belonging to nearby 
dwarf elliptical, dwarf spheroidal, spiral and giant elliptical 
galaxies (e.g., Buonanno et al. 1998; Layden & Sarajedini 
1997; Harris, Poole & Harris 1998; Cohen, Blakeslee & 
Rhyzov 1998; Puzia et al. 1999; Lee & Kim 2000; Beasley 
et al. 2000; Larsen et al. 2001). Indeed, evidence from 
various sources (see Gnedin et al. 2001) points to mini- 
mum GC ages in the range 8 ^ t <; 12 Gyr, suggesting 
that, in the currently fashionable ACDM models, the bulk 
of GC formation in early-type galaxies had occurred before 
1 ^ ^ ^ 4. Thus, for the early- type galaxies considered 
here, we may be confident that the observed GC colors are 
reliable metallicity indicators. 

Several recent studies have reported color/metallicity 
distributions for the GC systems of early-type galaxies 
based on archival HST data {e.g., Gebhardt & Kissler- 
Patig 1999; Kundu & Whitmore 2001; Larsen et al. 2001). 
As noted by van den Bergh (2001), the GC systems of 
these galaxies exhibit "a bewildering variety of character- 
istics" . This is particularly true of their metallicity distri- 
butions which range from obviously multi-modal systems 
to those with a single (metal-rich or metal-poor) peak. 
This remarkable variety has led some investigators to sug- 
gest that no single recipe for galaxy formation is capable 
of explaining the full range in GC properties {e.g., Har- 
ris 2000). Indeed, there is clear supporting evidence for 
most of the models that seek to explain the formation of 
GC systems and their host galaxies: e.g., the formation 
of massive clusters in gas-rich mergers {e.g., Whitmore & 
Schweizer 1995) and the acquisition of GCs through dissi- 
pationless accretions {e.g., Ibata, Gilmore & Irwin 1994). 
The most immediate question is which — if any — of these 
processes has dominated the formation of the GC systems 
which we observe today. 



In what follows, we model these diverse GC metallicity 
distributions using a simple prescription for hierarchical 
galaxy formation, described more fully in §3. In the in- 
terests of homogeneity, we limit ourselves to the GC color 
and metallicity distributions for 28 early-type galaxies pre- 
sented in Kundu & Whitmore (2001). Their galaxy sample 
is ~ 60% larger than that of Larsen et al. (2001), while 
the number of detected GCs per galaxy is ~ 50% larger 
than that of Gebhardt & Kissler-Patig (1999) for galaxies 
appearing in both studies. A single galaxy in the Kundu 
& Whitmore sample, NGC 4486B, has been excluded from 
the analysis since the modeling of its GC metallicity dis- 
tribution is complicated by the fact that a significant frac- 
tion of its initial mass has probably been tidally stripped 
by NGC 4486 {e.g., Faber 1973; Kormendy et al. 1997). 

Kundu & Whitmore (2001) present VI photometry and 
(V-I) colors for GC candidates associated with each of 
their program galaxies. The number of GCs, A^gc , ranges 
between 445 for NGC 4649 and 41 for NGC 5845, with a 
mean of (Nqc) = 163. Since our goal is to model the de- 
tailed shape of the observed GC metallicity distributions, 
we restrict ourselves to GCs with colors measured to a 
precision of cr(V-I) = 0.2 mag or better — roughly the 
separation between the metal-rich and metal-poor peaks 
in the most luminous galaxies. With this selection cri- 
terion, the number of GCs per galaxy ranges from 404 in 
NGC 4649 to 41 in NGC 5845, with a mean value of (A^gc) 
= 144. The adopted GC color-metallicity relation is that 
of Kissler-Patig et al. (1998), 

[Fe/H] = -4.50 + 3.27(V - I) (1) 

which is based on a combined sample of GCs in the Milky 
Way and NGC 1399 with measured metallicities. Al- 
though there is significant scatter about this best-fit re- 
lation (see Figure 13 of Kissler-Patig et al. 1998), our 
conclusions are largely insensitive to the exact choice of 
color-metallicity relation. The clusters used to establish 
this relation have metallicities in the range —2.4 ^ [Fe/H] 
^ -1-0.3 dex (0.64 ^ V-I ^ 1.47 mag), an interval which 
includes the vast majority of clusters considered here. 

To simulate the metallicity distribution of a GC system, 
we require an estimate of the host galaxy's luminosity, and 
hence distance. In principal, it is possible to measure in- 
dividual galaxy distances from the GC data themselves, 
using the observed turnover of the GC luminosity func- 
tion (GCLF). Unfortunately, Kundu & Whitmore (2001) 
were able to obtain reliable two-parameter GCLF fits {i.e., 
turnover magnitude and dispersion) for just 10 of their 28 
galaxies. As a consequence, distances for their program 
galaxies were taken primarily from early versions of the 
surface brightness fluctuation (SBF) survey of Tonry and 
collaborators, with extinction estimates from Burstein & 
Heiles (1984). Recently, the final version of the SBF cat- 
alog — based on an improved calibration — has become 
available (Tonry et al. 2001), so we adopt the latest SBF 
distances for the sample galaxies. Likewise, improved ex- 
tinctions for all 28 galaxies have been calculated from the 
Schlegel, Finkbeiner & Davis (1998) DIRBE maps. Un- 
fortunately, two of the sample galaxies — NGC 7626 and 
NGC 5982 — are not included in Tonry et al. (2001) cat- 
alog. For NGC 7626, we have combined the Schlegel et 
al. (1998) extinction with the Tully-Fisher distance for 
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the Pegasus group from Sakai et al. (2000) to find a true 
distance modulus of 33.73 mag. This is significantly larger 
than the value of 33.01 adopted by Kundu & Whitmore 
based on the fundamental plane distance of Prugniel & 
Simien (1996). For NGC 5982, no recent distance deter- 
mination is available, so we follow Kundu & Whitmore in 
adopting the Prugniel & Simien (1996) distance for this 
galaxy, but corrected for our new extinction. Because of 
these revised distances and extinctions, the absolute mag- 
nitudes of the program galaxies differ from those reported 
in Kundu & Whitmore. The mean magnitude difference is 
(AMy) = 0.18±0.08 mag (mean error), in the sense that 
our magnitudes are brighter than those of Kundu & Whit- 
more. The largest discrepancies are found for NGC 5813, 
NGC 7626 and IC 1459, whose absolute magnitudes have 
become brighter by AMy ~ 0.9 mag. 

Observational data for the sample galaxies are summa- 
rized in Table |l|. From left to right, the first six columns in 
this table report the galaxy name, integrated V-band mag- 
nitude from NED, extinction from Schlegel et al. (1998), 
adopted distance modulus, absolute visual magnitude, and 
total number of GCs having (t(V-I) < 0.2 mag. The re- 
maining columns are discussed in detail below. 

3. NUMERICAL METHOD 

3.1. Model Input 

We model the observed GC color and metallicity distri- 
butions for our program galaxies using a numerical code 
that has been described in two previous papers. The first 
of these papers (Cote et al. 1998) outlined the strategy 
for simulating GC color and metallicity distributions of 
galaxies that grow hierarchically via dissipationless merg- 
ers, and applied the method to the observed GC metallic- 
ity distribution of the giant elliptical galaxy NGC 4472, 
one of the objects in the current sample. A number of re- 
finements to this original code were described in Cote et al. 
(2000), where the model was extended to the GC systems 
and halo field stars of spiral galaxies {i.e., the Milky Way 
and M31). Since a complete discussion of the approach 
may be found in these papers, we give only a brief outline 
here. 

The method uses an empirical, Monte-Carlo approach 
to simulate the metallicity distributions of GCs associated 
with the end-products of multiple dissipationless mergers 
and accretions, as expected in hierarchical models in which 
the bulk of GC (and star) formation takes place in distinct 
protogalactic fragments and occurs prior to the epoch of 
galaxy assembly.^ Using a Monte-Carlo approach, it is pos- 
sible to simulate the GC metallicity distribution of an in- 
dividual galaxy by specifying three model ingredients: (1) 
an assumed form for the mass or luminosity function of 
the protogalactic fragments which agglomerate into the 
final galaxy; (2) a relation between the mass or luminos- 
ity of the protogalactic fragment and the mean metallicity 
of its GC system; and (3) the number of GCs per unit 
protogalactic fragment mass or luminosity. Each of these 
ingredients is determined empirically from observations of 
the GC systems of low-redshift galaxies. This empirical 
approach has the important practical advantage that it 
avoids the uncertainties involved in modeling directly the 

^Strictly speaking, the process is dissipationless only in the sense thi 
by the mergers themselves, but rather, by the cooling and collapse of g 



cycle of gas cooling, star formation and feedback in these 
halos (issues which, ultimately, must be addressed in any 
complete theory of GC formation). As discussed in Cote 
et al. (1998), the merging probability of protogalactic frag- 
ments is assumed to be independent of mass. As Figure 1 
of Bullock et al. (2001) shows for the specific case of the 
Milky Way, this assumption appears remarkably sound for 
protogalactic fragments with circular velocities of less than 
Vc ~ 60 km s~^, corresponding to a mass of ~ a few 
times 10^ Mq- 

The luminosity function of protogalactic fragments is 
parameterized as a Schechter function, 

n{C) oc (£/£*)" exp {-C/C*) (2) 

where C* is a characteristic luminosity and a is an ex- 
ponent which governs the relative number of faint and 
bright systems (Schechter 1976). Note that the luminos- 
ity in equation (2) refers to the present day luminosity 
of surviving protogalactic fragments (i.e., isolated galax- 
ies of low- and intermediate- luminosity; Larson 1988). 
Over the luminosity range considered here, we assume a 
constant F-band mass-to-light ratio of Ty = 5 at z = 0, 
which is roughly appropriate for the low- and intermediate- 
luminosity galaxies used below to define the relation be- 
tween mean GC metallicity and host galaxy luminosity 
(Mateo 1998). We take the low-luminosity cutoff in equa- 
tion (2) to be Cv = 1.6 x 10^ Cv,e, or 7W = 8 x 10^ Mq. 
These limits correspond to the Fornax dSph galaxy (the 
faintest galaxy known to contain its own GC system) but 
the precise choice has little effect on the calculated GC 
metallicity distributions for slopes shallower than a — —2. 
At the bright end, we adopt a constant value of jCy = 
2 X 10^° ^v,Q for all simulations (Cote et al. 2000), while a 
is a free parameter to be determined on a galaxy-by-galaxy 
basis through a comparison of the observed and simulated 
GC metallicity distributions. Note that, for masses less 
than A4 ^ 10^^ -Mq, the protogalactic mass function has 
a nearly power-law form, 

n{M) oc {M/M*r, (3) 

which is also characteristic of the extended Press-Schechter 
mass functions used in semi-analytic models. Our decision 
to treat the change in the initial luminosity and mass func- 
tions as a change in a alone is clearly an approximation, as 
C* and A4* must have unavoidably increased in the pres- 
ence of mergers and accretions. Although our simulations 
have no time resolution, cosmological N-body simulations 
demonstrate that the change in 7W* should have only a 
modest effect on the final GC metallicity distribution: Fig- 
ure 12 of Pearce et al. (2001) suggests the that change in 
the high-end cutoff of the mass function between z = 3 
and z = 0.2 is about a factor of four, corresponding to a 
maximum shift in (V-I) color of 0.08 mag for the highest 
mass protogalactic fragments. 

The relation between GC metallicity and luminos- 
ity/mass of the host protogalactic fragment has been dis- 
cussed extensively in Cote et al. (2000). Figure |^ shows 
the mean GC metallicity plotted against the luminosity (or 
mass) of nine dwarf elliptical (dE) and dwarf spheroidal 

it the formation of GCs — clearly a dissipative process — is not driven 
as with each protogalactic fragment. 



4 



Globular Clusters and the Missing Satellite Problem 



(dSph) galaxies (circles), with data taken directly from 
that paper. Note that GC metallicities for these galaxies 
have been converted into (V-I) colors according to equa- 
tion (1). Triangles indicate the metal-rich GC systems of 
two well-studied spiral galaxies: the Milky Way and M31. 
In both cases, the luminosities and masses refer to their 
bulges; that is to say, we assign the metal-rich GCs in these 
galaxies to their bulge components, and not their disks 
(Minniti 1995; Cote 1999). Our estimate for the relation 
between the mean GC color and the absolute magnitude 
of the host protogalactic fragment is 

(V - I) = 0.04 - 0.055My, (4) 

which is indicated by the dashed line in Figure |l|. This 
"zero-age" relation gives the dependence of mean GC color 
on the luminosity, or mass, of the progenitor protogalac- 
tic fragment, as it would appear today, if the subsequent 
evolution has been mostly passive. In the simulations, the 
mean color/metallicity of GCs belonging to each proto- 
galactic fragment is determined using equation (4), with 
an assumed la scatter of 0.03 mag (0.1 dex). The inter- 
nal dispersion in color for each protogalactic fragment is 
also approximated as a Gaussian, with dispersion 0.09 mag 
(0.29 dex) based on the findings of Cote et al. (2000). 

Finally, we adopt a constant GC specific frequency for 
all protogalactic fragments — Sn = 5 ± 1 at the present 
time — and neglect possible variations in Sn at the low- 
and high-mass ends of the galaxy distribution (see §6). In 
the notation of McLaughlin (1999), this specific frequency 
corresponds to a GC formation efficiency of 

ecc = 0.29 ± 0.06% (5) 

by mass. This is consistent with the value of ecc = 
0.26±0.05% advocated by McLaughhn (1999) based on ob- 
servations of about one hundred early-type galaxies. The 
20% uncertainty in our adopted values of Sn and ed re- 
flects the observed scatter about their means, and has been 
included explicitly in our Monte-Carlo simulations. 

3.2. An Example 

Before proceeding, we present a descriptive example of 
how the GC system of a galaxy that grows hierarchi- 
cally through dissipationless mergers will evolve in the 
color-magnitude plane. Consider a protogalactic fragment, 
A, with absolute magnitude My = —20 mag and mass 
Ai — 4: X 10^" Mq which falls on the zero-age color- 
magnitude relation (i.e., the dashed line in Figure |l|). Ac- 
cording to equation (4), the unimodal GC system associ- 
ated with this protogalactic fragment will have (V-I) ~ 
1.14 mag. First, suppose that this fragment merges with 
an identical fragment, producing an end-product, CI, with 
My = 20.75 mag and a GC system having a single peak 
at (V-I) = 1.14 mag. Alternatively, fragment A could ac- 
crete a population of 100 protogalactic fragments, B, each 
of which has My = — 15 mag, and (V-I) = 0.865 mag. But, 
by virtue of equation (5) and our assumed constant mass- 
to-light ratio, fragment A contributes the same number of 
GCs as the 100 smaller fragments combined. Thus, the 
GC systems of the final galaxy, C2, will have two distinct 
components (containing equal numbers of metal-rich and 
metal-poor GCs), and a mean color of (V-I) ~ 1.00 mag. A 



third possible end-product, C3, is also shown in Figure ^ 
In this case, the galaxy was formed entirely through the 
agglomeration of (200) fragments of type B, so that the 
final GC system is unimodal, but with a metal-poor peak 
at (V-I) = 0.865 mag. 

The point of this simple exercise to emphasize that, once 
the "zero-age" color-metallicity is fixed, it is possible — 
through a suitable adjustment of the mass spectrum of the 
protogalactic fragments — to reproduce the GC systems 
of any galaxy whose GC system falls within the dotted 
region shown in Figure |l|. This simple recipe for galaxy 
formation combines dissipationless mergers, which serve to 
move the galaxy to the right of the initial color-magnitude 
relation, with a dissipative, monolithic collapse scenario 
for the constituent protogalactic fragments. Our goal is to 
use the observed GC metallicity distribution for a specific 
galaxy and solve an inverse problem to determine the mass 
spectrum of protogalactic fragments which best matches 
the observed metallicity distribution. 

3.3. Application 

With equations (2) and (4) specified, it is a trivial mat- 
ter to calculate analytically the metallicity distribution of 
a GC system which arises from dissipationless mergers and 
accretions. In practice, such an approach is thwarted by 
the highly stochastic nature of hierarchical galaxy forma- 
tion {e.g., Lacey & Cole 1993), particularly at the high- 
mass end of the protogalactic mass spectrum. Fortunately, 
this stochasticity can be used to reconstruct merger his- 
tories for individual galaxies since the fossil record of the 
largest protogalactic fragments is ussually preserved in the 
final GC metallicity distribution. 

To do so, we generate, for each galaxy in our sample, 50 
simulated GC systems at each point over a 10x10 grid in 
(C, ct) where C is the ratio between the mass of the largest 
protogalactic fragment and the final mass of the galaxy, 

C = M'/Mf, (5) 

and a is the slope of the protogalactic mass function. The 
(C, a) plane is sampled at ^ = 0.1, 0.2, ...1.0 and and 
a — —1.1, —1.2, ... — 2.0. For each of the 5000 simulations 
per galaxy, we extract at random from the full GC system 
the same number of GCs as contained in the Kundu & 
Whitmore (2001) catalog, add the appropriate amount of 
measurement error to the simulated (V-I) colors, and bin 
the actual and simulated data in an identical manner (see 
§4.1). A goodness-of-fit statistic, 

is then calculated for each simulation, along with a mean 
for each point on the 10x10 grid. The number of bins, 
TVbin, is chosen on a galaxy- by-galaxy basis as described 
in §4.1. As an illustration of the method, Figure || shows 
a contour plot over this grid for one the galaxies in our 
sample, NGC 4472. From these simulations, we deduce 
best-fit values of C = 0.30±0.05 and a = -1.85±0.10 {la 
uncertainties) . The cross in Figure || shows the location 
in the (C, a) plane of the simulation which produced the 
best match to the observed GC metallicity distribution. 
These results are summarized in the final three columns 
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of Table |i| which record the for the best-fit simulation, 
along with the corresponding values of a and C and their 
1(7 uncertainties. 

The mean simulated color of the GC systems in the 
sample is (V-I)sim = 1.00, while the mean difference be- 
tween the observed and simulated GC color is (V-I)obs-(V- 
I)sim — 0.00 mag, with a standard deviation of 0.01 mag. 
We now turn our attention to a comparison of the detailed 
shapes of the observed and simulated GC color and metal- 
licity distributions. 

4. RESULTS 

4.1. Metallicity Distributions 

In this section, we compare the observed GC color and 
metallicity distributions with the best-fit models for each 
of our program galaxies. We follow the usual practice of 
presenting the results as histograms: i.e., as parametric 
representations of the underlying distributions. An opti- 
mal histogram bin width, h, is calculated on a case-by-case 
basis according to the GC sample size and level of non- 
Gaussianity for each metallicity distribution. Specifically, 
we use the method of Scott (1992) and set the bin width 
at 



N, 



1/3 



(6) 



GC 



where a is the standard deviation of the distribution and 
77 is a factor in the range 0.4 <C, 77 <^ 1 related to the dis- 
tribution's skewness. For our program objects, the bin 
widths range between h ~ 0.055 mag for NGC 4649, and 
h ~ 0.138 for NGC 4550; these values of h are comparable 
to the mean errors, a{V-T) = 0.08 mag, of the measured 
GC colors. Note that this procedure differs from that of 
Kundu & Whitmore (2001), who opted to fix the bin width 
at a constant value of ft. = 0.03 mag, although A^gc varies 
by more than an order of magnitude among their program 
galaxies. 

Figures ||-^ show the metallicity distributions for GCs 
belonging to the 28 galaxies listed in Table |^. The ob- 
served distributions are indicated by the filled circles, while 
the errorbars show the Poisson uncertainties in each bin. 
For each galaxy, the best-fit model is indicated by the 
dashed curve, which has been binned in the same man- 
ner as the actual data. The excellent agreement between 
the observed and simulated distributions is noteworthy in 
light of the fact that, for the simple model considered here, 
the effects of just two free parameters, C. and a, are ex- 
plored. A more comprehensive exploration of parameter 
space might include possible mass-dependent variations in 
the GCs formation efficiency, the dynamical evolution of 
the GC systems, GC formation induced by the merger and 
accretion process, and the form of the intrinsic metallic- 
ity distribution in the protogalactic fragments. The first 
two processes were included in the Milky Way formation 
model discussed in Cote et al. (2000), but have been omit- 
ted here for the sake of simplicity. We note, however, that 
both the weak dependence of GC specific frequency on host 
galaxy luminosity {e.g., Durrell et al. 1996; Miller et al. 
1998; McLaughhn 1999) and the preferential disruption of 
the more centrally concentrated metal-rich GCs in these 



galaxies, will serve to increase the ratio of metal-poor to 
metal-rich GCs. Including these effects in the simulations 
would give additional freedom to match the relative num- 
bers of metal-rich and metal-poor GCs in these galaxies. 
To first order, these processes affect the peak heights for 
the subsystems, and not their mean metallicities. 

Figure ^ presents the same color-magnitude plane for 
GC systems that was shown in Figure H. However, we 
now include the GC systems of the 28 galaxies listed in 
Table |]. Kundu and Whitmore (2001) classified each of 
these GC systems as either "certainly" bimodal, "likely" 
bimodal or unimodal. We adopt their classifications ver- 
batim. For the 16 galaxies classified as certainly or likely 
bimodal, the peak color for the metal-rich and metal-poor 
GC systems (according to KMM) are indicated by the 
filled squares and circles, respectively.^ The filled triangles 
indicate the mean color and metallicity of the 12 unimodal 
GC systems. The open symbols show the median colors 
and metallicities of the GC systems in each of these galax- 
ies based on 100 simulations using the best-fit values of ( 
and a listed in Table ^ Simulations for the unimodal sys- 
tems are indicated by the open triangles, while the open 
squares and open circles show the median colors for the 
metal-rich and metal-poor components in the bimodal sys- 
tems (determined from the best- fitting double Gaussians). 
To illustrate the typical range in simulated GC color and 
metallicity, an errorbar is included in the lower right cor- 
ner of this figure, showing the la dispersion in the colors 
and metallicities based on the simulations. 

The simulations reproduce the observed GC colors and 
metallicities of not just the unimodal systems, but also 
of the metal-rich and metal-poor components in the bi- 
modal GC systems. Likewise, the simulations reproduce 
the observed trend for the metal-poor peak in the bimodal 
GC systems to depend only weakly on host galaxy lumi- 
nosity (Forbes, Brodie & Grillmair 1997; Cote et al. 1998; 
Larsen et al. 2001) and, by definition, the clear correlation 
between metal-rich peak and host galaxy luminosity. This 
latter trend can be traced to the fact that the first-ranked 
protogalactic fragments: (1) themselves originate on the 
"zero-age" color-magnitude relation; and (2) contribute a 
significant fraction of the galaxy's final mass. Indeed, for 
each galaxy classified by Kundu & Whitmore as bimodal, 
we show the initial location in the color-magnitude plane 
of the largest protogalactic fragment, determined from the 
best-fit simulation for each galaxy (open stars). Thin 
arrows connect these initial positions to the locations of 
the metal-rich peak measured directly from the observed 
GC metallicity distribution. Note that, in a few cases, 
measurement errors and observational scatter (indicated 
by the errorbar) cause the initial location in the color- 
magnitude plane to lie below the final location of the metal- 
rich peak. The filled pentagons show the GC systems of 
the same protogalactic fragments (i.e., dwarfs and spiral 
bulges) used to define the "zero-age" color-magnitude re- 
lation in Figure |[ reproduced as the dashed line in Fig- 
ure |l^. We conclude that the simple model considered 
here provides an adequate match to both the observed 
GC metallicity distributions and the final location of the 
GC systems in the color-magnitude plane. 



^The GC systems of NGC 3377 and NGC 584 are classified by Kundu & Whitmore (2001) as "likely" to be bimodal in their Table Q 
but plotted as "certain" in their Figure 5. We have assumed that the GC systems of both galaxies show statistically significant evidence of 
bimodality. 
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4.2. Comparisons with Previous Results 

On the basis of the same GC data analyzed here, Kundu 
& Whitmore (2001) claim that hierarchical models for 
galaxy formation are unviable. As Figures p[-[To| demon- 
strate, this conclusion is not supported by their data. 

Kundu & Whitmore (2001) made no attempt to sim- 
ulate the GC color distributions expected in hierarchical 
formation scenarios, despite the fact that the implementa- 
tion of such algorithms has been described in detail previ- 
ously {e.g., Cote et al. 2000; Harris et al. 2000). Instead, 
they advocate a qualitative scenario in which GC forma- 
tion occurs in multiple star formation episodes driven by 
mergers of gas-rich disk galaxies, although such models 
are curently not amenable to observational tests. They 
further note that some of the low-luminosity galaxies in 
their sample {e.g., NGC 1439, NGC 3377 and NGC 4660) 
are likely to have bimodal GC metallicity distributions and 
claim that this observation poses "a serious problem" for 
the hierarchical model. This claim is surprising since Cote 
et al. (2001) showed that hierarchical models were able 
to reproduce the bimodal GC metallicity distribution of 
the Galactic spheroid, which is several times fainter than 
NGC 1439, NGC 1427 and NGC 3377. Moreover, all of 
the low-luminosity galaxies having bimodal GC metallic- 
ity distributions fall within the region in Figure |^ that is 
expected to be occupied by galaxies which form hierarchi- 
cally (§3.2). As Figures attest, the observed diversity 
in the GC metallicity distributions and the bimodal na- 
ture of some low-luminosity ellipticals could arise naturally 
within the context of a simple hierarchical model. 

Kundu & Whitmore (2001) also state that the success 
of the hierarchical model "is critically dependent on the 
second-order nature of [the metallicity-luminosity relation 
for] the systems being cannibalized" . However, this claim 
too is incorrect since the model distributions presented in 
Figures 0-0 are derived from a linear metallicity-luminosity 
relation (see §3.1), not the second-order relation adopted 
by Cote et al. (1998) on these basis of preliminary data. 
While it is clearly important to further refine the metal- 
licity luminosity relation for the GC systems of low- and 
intermediate-luminosity galaxies, the model described here 
is quite general in nature, and our basic conclusions are 
unlikely to be changed by minor modifications to the 
metallicity-luminosity relation. 

4.3. Protogalactic Mass Spectra 

Having shown that hierarchical growth through dissipa- 
tionless mergers and accretions offers a viable explanation 
for the diverse metallicity distributions of GC systems in 
early-type galaxies, we now examine properties of proto- 
galactic fragments from which the galaxies may have been 
assembled. Figures |l^0 show protogalactic luminosity 
and mass functions for each galaxy, based on the simula- 
tion which matches the observed GC metallicity distribu- 
tion. The shaded histogram in each panel shows the best 
estimate for the distribution of protogalactic fragments, 
while the dashed curve indicates the parent Schechter func- 
tion from which the simulation was drawn. For all galax- 
ies, the corresponding GC metallicity distributions are in- 
dicated by the dashed curves in Figures ^]9[ 

A striking feature of the protogalactic luminosity and 
mass functions shown in Figures [TlfjlTl is their steep 



slope. From Table yj, the average power-law exponent 
for the protogalactic luminosity and mass functions is 
(a) — —1.88 ± 0.03 (mean error), with a standard devi- 
ation of cr(a) — 0.13. This value is consistent with previ- 
ous findings for NGC 4472 {a 1.8; Cote et al. 1998), 

the Milky Way (a 2.0; Cote et al. 2000) and M31 

{a ~ —1.8; Cote et al. 2000). In the context of hier- 
archical formation models, it is the apparent universality 
of this exponent which, when coupled with the fixed low- 
mass cutoff of the protogalactic luminosity function, con- 
spires to produce the population of metal-poor GCs with 
[Fe/H] —1.2 dex that are present in virtually every large 
galaxy studied to date {e.g.. Ashman & Bird 1993; Brodie, 
Larsen & Kissler-Patig 2000). As the same time, the ob- 
served diversity of GC metallicity distributions arises from 
the stochastic fluctuations about this "universal" proto- 
galactic mass spectrum, particularly at the high-mass end. 
By virtue of their large mass, these fragments contribute 
the bulk of the metal-rich GCs associated with the final 
galaxy. 

From Table ||, we see that there is a fairly wide range 
in the properties of the most massive protogalactic frag- 
ments. On average, the most massive fragment contributes 
a fraction {() = 0.25 ± 0.03 (mean error) of the galaxy's 
final mass, with a standard deviation of a{Q — 0.14. For 
some galaxies, such as NGC 4406, the contribution of the 
largest fragment amounts to no more ~ 10% of the final 
mass, while in NGC 4458, as much as 70% of the final mass 
may have arisen in the dissipational collapse of the largest 
protogalactic fragment. For NGC 4472, we find a best-fit 
value of C = 0.30 ± 0.05 (mean error), which compares 
favorably with the value of C ~ 0.35 estimated by Cote 
et al. (1998) using ground-based observations of its GC 
system carried out in a different filter system {i.e., CTi 
photometry from Geisler, Lee & Kim 1996), and using an 
alternative color-magnitude relation. 

Recent studies have concluded that approximately 50% 
of bright, early- type galaxies show evidence for bimodal 
color distributions (Gebhardt & Kissler-Patig 1999; Larsen 
et al. 2001; Kundu & Whitmore 2001), so it is of inter- 
est to ask how often the simulations yield a metallicity 
distribution that is bimodal at a statistically significant 
level. Unfortunately, the answer to this question depends 
rather sensitively on the properties of each dataset used 
to assess the level of biniodality. To understand why this 
is so, consider Figure |l^, which shows the metallicity dis- 
tribution for a representative program galaxy, NGC 5322. 
For the purposes of illustration, the solid curve shows the 
(normalized) metallicity distribution corresponding to the 
protogalactic mass function plotted in the upper panel of 
Figure |l^, where the intrinsic dispersion in metallicity for 
each protogalactic fragment is taken to be 0.01 dex. As dis- 
cussed in §3.1 and in Cote et al. (2000), actual fragments 
are likely to have had intrinsic dispersions much larger 
than this, but this extreme example serves to illustrate 
the complex nature of the underlying metallicity distribu- 
tion, and demonstrates the dangers of interpreting multi- 
modal metallicity distributions as prima facie evidence for 
multiple bursts of GC formation. The remaining curves 
show this same (renormalized) distribution after smooth- 
ing with Gaussians having dispersions of 0.1, 0.2, 0.3 and 
0.4 dex (intended to represent the combined effects of an 
intrinsic metallicity dispersion and observational errors). 
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This last distribution is, like the observed distribution, os- 
tensibly unimodal in nature (Kundu & Whitmore 2001), 
but this is mainly a consequence of the intrinsic nietallic- 
ity dispersion (~ 0.3 dex) for the protogalactic fragments, 
and the ^ 0.3 dex dispersion associated with metallicity 
measurements derived from (V-I) colors having a preci- 
sion of roughly 0.1 mag. Thus, establishing the frequency 
of bimodality depends rather sensitively on the precision 
of metallicity measurements (and on other factors includ- 
ing the total luminosity of the final galaxy and its detailed 
merger tree). 

5. IMPLICATIONS FOR THE MISSING SATELLITE 
PROBLEM 

As mentioned in §1, one of the most troubling prob- 
lems facing CDM models on small scales is the so-called 
"missing satellite" problem. Simply stated, the number 
of low-mass dark halos predicted by CDM models exceeds 
the observed number of dwarf galaxies by factors of tens, 
or even hundreds {e.g., Kauffmann, White & Guiderdoni 
et al. 1993; Klypin et al. 1999; Moore et al. 1999). Is it 
possible to explain the underabundance of dwarf galaxies 
in the local universe, and at the same time retain the basic 
framework of CDM models? 

Figure |l^ shows the ensemble spectrum of protogalactic 
fragments in our sample: i.e., summed over all 28 galaxies 
(filled circles). As with the individual galaxies, the com- 
posite mass function is found to follow a power-law behav- 
ior, n{Ai) oc Ai°', with exponent a ~ —2. While this mass 
function is much steeper than the local galaxy luminosity 
function {i.e., Pritchet & van den Bergh 1999; Folkes et al. 
1999; Blanton et al. 2001), it is nevertheless in remark- 
able agreement with the mass spectrum of dark matter 
halos predicted by CDM models of structure formation; 
the dotted line in Figure [l^ shows the analytical approx- 
imation, n{A4) oc M.~^, for the mass spectrum of dark 
matter halos in the N-body simulations of Moore et al. 
(1999). In short, we find evidence for a disrupted popula- 
tion of dwarf-sized protogalactic fragments in each of our 
program objects, whose numbers and masses closely resem- 
ble those of the predicted population of ''missing" satellite 
galaxies. 

This finding may have important consequences for the 
various models which have been advanced as possible so- 
lutions to the missing satellite problem. Such a disrupted 
population is, in fact, predicted by models which account 
for the missing satellites by suppressed gas accretion in 
low-mass halos after the epoch of reionization due the 
presence of a strong photo-ionizing background {e.g., Rees 
1986; Thoul & Weinberg 1996; Bullock, Kravtsov & Wein- 
berg 2000). A disrupted component of this sort is also ex- 
pected in models which seek to explain the missing satel- 
lite problem through self- interacting CDM {e.g., Spergel 
& Steinhardt 2000; Dave et al. 2001). By contrast, a dis- 
rupted component is not expected in models which explain 
the underabundance of dwarf galaxies through an ah ini- 
tio reduction in small scale power {e.g., Kamionkowski & 
Liddle 2000) or by suppressing the formation of low-mass 
halos via warm dark matter {e.g., Colin, Avila- Reese & 
Valenzuela 2000; Sommer-Larsen & Dolgov 2001). Thus, 
these models appear to be inconsistent with the observed 
metallicity distributions of GCs in galaxies at z ~ since 
the same low-mass halos whose formation is suppressed 



in these scenarios are actually required to explain the GC 
metallicity distributions of luminous galaxies in hierarchi- 
cal cosmologies. 

6. SUMMARY AND CONCLUSIONS 

We have described a technique for simulating the metal- 
licity distributions of GCs belonging to galaxies that grow 
hierarchically through dissipationless mergers and accre- 
tions. Applying the model to the GC systems of 28 early- 
type galaxies for which high quality HST data are avail- 
able demonstrates such a mechanism is capable of explain- 
ing quantitatively the observed diversity in extragalactic 
GC systems. Thus, claims that the metallicity distribu- 
tions of GCs in early-type galaxies are inconsistent with a 
hierarchical origin (Kundu & Whitmore 2001) are unsub- 
stantiated. 

Protogalactic mass spectra for the sample of early-type 
galaxies considered here reveal a striking similarity: they 
appear to follow a power-law distribution in mass, with 
n{M.) cx A4~^. This mass spectrum is indistinguishable 
from the mass spectrum of dark matter halos predicted by 
CDM models, and thus has an important implication for 
the missing satellite problem (Klypin et al. 1999; Moore 
et al. 1999). Specifically, we find evidence from the GC 
systems of the galaxies examined here for a disrupted pop- 
ulation of dwarf-sized protogalactic fragments whose num- 
bers and masses closely resemble those of the predicted 
population of missing satellite galaxies. 

If the basic scenario outlined here is correct {i.e., the 
early formation of GCs in low- and intermediate-mass ha- 
los after a rapid period of mass-dependent chemical enrich- 
ment, followed by the agglomeration of these halos into 
the main body of the galaxy), then there may be some 
interesting observational consequences of this picture. For 
instance, if most of the low-mass dark matter halos pre- 
dicted by CDM models have indeed been rendered unde- 
tectable as a result of disruption and the suppression of 
gas accretion after reionization, as suggested by Bullock 
et al. (2000), then some fraction of the original halo pop- 
ulation may survive as "dark galaxies" (see, e.g., Klypin 
et al. 1999 and references therein). Precise estimates for 
the redshift of GC formation are not yet available, but 
current evidence seems to favor formation at 2: ^ 7 (see 
Gnedin et al. 2001), close to the expected edge of the 
epoch of reionization (Djorgovski et al. 2001; Becker et al. 
2001; Barkana 2001). This raises the possibility that some 
of these putative "dark galaxies" , which have experienced 
little or no star formation, might nevertheless have man- 
aged to form GCs. Interestingly, there is some evidence 
that the number of GCs per unit host galaxy luminos- 
ity increases toward fainter systems, at least among dE,N 
galaxies (Durrell et al. 1996; Miller et al. 1998). Searches 
of extragalactic GC systems for compact, dynamically- 
bound groupings of GCs — presumably metal-poor in na- 
ture and not associated with an identifiable satellite galaxy 
— might prove interesting. 

We thank Arunav Kundu for providing the GC pho- 
tometry used in the analysis, and an anonymous ref- 
eree for helpful comments. MJW was supported by NSF 
grant ASTOO-71149. This research has made use of the 
NASA/IPAC Extragalactic Database (NED) which is op- 
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Table 1 

Protogalactic Mass Functions of Early-Type Galaxies Deduced From Their Globular Cluster 

Systems 
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.27 


185 


0, 


.96 


^2, 


.00 ±0.05 


0, 


.10 ±0.05 


NGC 4589 


10 


.73 


0, 


.038 
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.80 


-21, 


.07 
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1, 


.03 
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0, 
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0, 


.038 
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.07 
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.87 
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0, 


.44 


-2, 


nn-HO.lO 
■""-0.05 


0, 


.30 ±0.10 


NGC 3379 


9, 


.28 


0, 


.033 


31, 


.20 


-20, 


.92 


63 


0, 


.44 


-1, 


qn-fO.15 
■^"-0.10 


0, 


.35 ±0.10 


NGC 821 


10 


.68 


0, 


.143 


32, 


.26 


-21, 


.58 


101 


1, 


.19 


-2, 


nn+O.lO 
■""-0.05 


0, 


in+0.10 
■-•■"-0.05 


NGC 3608 


10 


.76 


0, 


.029 


31, 


.87 


-21, 


.11 


95 


1, 


.29 


-1, 


OC,+0.20 
■^^-0.10 


0, 


.10 ±0.05 


NGC 4291 


11 


.47 


0, 


.050 


32 


.21 


-20, 


.74 


138 


0, 


.84 


-1, 


qr+0.15 

■^^-0.10 


0, 


.10 ±0.05 


NGC 1439 


11 


.39 


0, 


.041 


32, 


.23 


-20, 


.84 


85 


1, 


.08 


-1, 


■^U-0.20 


0, 


lA + O.lO 

■-■^"-0.05 


NGC 1427 


10 


.86 


0, 


.018 


31, 


.90 


-21, 


.04 


146 


0, 


.54 


-2, 


QQ + 0.20 

■""-0.05 


0, 


9n+o.io 

■^"-0.15 


NGC 3377 


10 


.38 


0, 


.046 


30 


.36 


-19, 


.98 
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0, 


.35 


-1, 


on+0.45 
■^U_o.i5 


0, 
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NGC 4550 


11 


.68 


0, 


.057 


31, 


.13 


-19, 


.45 


47 


0, 


.83 


-1, 


7C-I-0.50 
■'^-0.30 


0, 


.20 ±0.15 


NGC 5845 


12 


.48 


0, 


.073 


32 


.24 


-19, 


.76 


41 


0, 


.35 


-1, 


7C-I-0.35 
■'^-0.30 


0, 


or-l-0.25 
■"5^-0.15 


NGC 4660 


11 


.24 


0, 


.048 


30 


.64 


-19, 


.40 


98 


0, 


.53 


-1, 


.65 ±0.15 


0, 


oc-l-O.lO 
■^3-0.05 


NGC 4458 


12 


.07 


0, 


.033 


31, 


.26 


-19, 


.19 


43 


0, 


.44 


-1, 


.75 ±0.25 


0, 


70+0.10 
■'^-0.05 



Mean = -1.88 ± 0.03 0.25 ± 0.03 



'^TuUy Fisher distance for the Pegasus group from Sakai et aL (2000), ApJ, 529, 698. 
'^Fundamental Plane distance from Prugniel & Simien (1996), A&A, 309, 749 
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Fig. 1. — Mean color and metallicity of globular clusters belonging to low- and intermediate-luminosity galaxies, plotted as a function of 
absolute galaxy magnitude (circles). Triangles show the metal-rich globular cluster populations of the Milky Way and M31, plotted against 
their respective bulge luminosities. The tightness of this correlation indicates that the chemical enrichment of the globular clusters in these 
systems was controlled primarily by the depth of the gravitational potential well in which they formed. The dashed line indicates our 
adopted color-magnitude relation for the GC systems of these protogalactic building blocks, each of which is assumed to have formed via a 
"monolithic" collapse. The dotted region shows the area occupied by galaxies which grow hierarchically via dissipationless agglomeration of 
these protogalactic fragments. 
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Fig. 2. — Confidence intervals for a and C> based on 5,000 Monte Carlo simulations of the metallicity distribution of globular clusters 
associated with the elliptical galaxy NGC 4472. In this notation, the slope of the protogalactic mass/luminosity function is denoted by a, 
while 1^ is the ratio between the mass of the largest protogalactic fragment and the final mass of the galaxy. The contours denote 68%, 90%, 
95% and 99% confidence limits. The cross indicates the location in the a-C plane of the simulation which produced the best match to the 
observed distribution. 



12 Globular Clusters and the Missing Satellite Problem 



[Fe/H] 

-2 -1 



60 



O 40 



CD 

^ 20 



"1 — I — I — r 



N4472 



4 



•4 



i 1 



_l I I I I \ I I I I I I I I I I L 



0.6 0.8 1 1.2 1.4 

(v-i) 

[Fe/H] 
-2 -1 



40 



m 
O 
O 



30 



I 20 

S 

^ 10 



T 1 1 1 1- 

- N4406 



1 — I — I — I — I — I — I — I — r 
4 



J I I I I I I L_ 



1 ^ 



0.6 0.8 1 1.2 1.4 

(V-I) 



[Fe/H] 

-2 -1 



60 



!/3 

u 



^ 40 
u 

CD 

B 

:i 20 



"T — I — I — I — r 



N4649 



i 



i 



■ •1 



^-•|- I I I 

_l I \ I I I 1 I \ 1 I \ I I I I L_ 



0.6 0.8 1 1.2 1.4 
(V-I) 



[Fe/H] 
-2 -1 







40 - 



5h 

B 20 
2; 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 
N4365 



1 



_l I I L 



0.6 0.8 



i-i 



1.2 1.4 



(V-I) 



Fig. 3. — Color anci metallicity distributions of globular clusters belonging to the early-type galaxies NGC 4472, NGC 4649, NGC 4406 
and NGC 4365 (filled circles). The errorbars show the Poisson uncertainties in each bin. The dashed curve shows the best- fit model for each 
galaxy, determined from 5,000 Monte Carlo simulations spanning a 10x10 grid in (a, Q. 
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Fig. 4.— Same as in Figure |, except for NGC 5322, NGC 4494, NGC 7626 and NGC 5982. 
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Fig. 5.— Same as in Figure |, except for IC 1459, NGC 3610, NGC 584 and NGC 4621. 
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Fig. 6.— Same as in Figure |, except for NGC 4552, NGC 5813, NGC 4589 and NGC 4278. 
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Fig. 7.— Same as in Figure |, except for NGC 4473, NGC 3379, NGC 821 and NGC 3608. 
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Fig. 8.— Same as in Figure |, except for NGC 4291, NGC 1439, NGC 1427 and NGC 3377. 
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Fig. 9.— Same as in Figure |, except for NGC 4458, NGC 5845, NGC 4660 and NGC 4458. 
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Fig. 10. — Mean color (mctallicity) of globular cluster systems, plotted against galaxy magnitude (mass). The dashed line shows the 
"zero-age" relation between galaxy mass and globular cluster mctallicity shown in Figure hi for unmerged protogalactic fragments: i.e., 
isolated galaxies of low and intermediate mass, and the bulge components of spiral galaxies (mled pentagons). The filled squares and circles 
show the observed colors for the metal-rich and metal-poor globular clusters in the 16 galaxies classified as bimodal by Kundu & Whitmore 
(2001). Open squares and circles indicate the median color of these components based on 100 simulations of the globular cluster metallicity 
distributions using the (a, C) values reported in Table hi. The filled and open triangles show the observed and simulated colors of globular 
clusters belonging to the 12 galaxies classified by Kundu & Whitmore (2001) as unimodal. The errorbar in the lower right corner shows the 
Icr scatter about these median values. For each galaxy with a bimodal metallicity distribution, a thin arrow connects the original position of 
the most massive protogalactic fragment (open stars) to the measured position of the metal-rich peak. 
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Fig. 11.— Mass and luminosity functions of protogalactic fragments for NGC 4472, NGC 4649, NGC 4406 and NGC 4365. The histograms 
show the distribution of protogalactic fragments for the best-fit models shown in Figure H. The dashed line in each panel indicates the parent 
Schechter function from which the protogalactic fragments were drawn. 
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Fig. 12.— Same as in Figure except for NGC 5322, NGC 4494, NGC 7626 and NGC 5982. 
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Fig. 13.— Same as in Figure except for IC 1459, NGC 3610, NGC 584 and NGC 4621. 
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Fig. 14.— Same as in Figure except for NGC 4552, NGC 5813, NGC 4589 and NGC 4278. 
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Fig. 15.— Same as in Figure except for NGC 4473, NGC 3379, NGC 821 and NGC 3608. 
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Fig. 16.— Same as in Figure except for NGC 4291, NGC 1439, NGC 1427 and NGC 3377. 
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Fig. 17.— Same as in Figure [l|, except for NGC 4458, NGC 5845, NGC 4660 and NGC 4458. 
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Fig. 18. — Simulated metallicity distribution for a representative galaxy, NGC 5322, assuming the intrinsic metallicity dispersion of each 
protogalactic fragment is 0.01 dex (solid curve). The remaining curves show this same distribution after smoothing with Gaussians having 
dispersions of 0.1, 0.2, 0.3 and 0.4 dex. 
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Fig. 19. — Comparison of the halo mass function measured from cosmological N-body simulations (dotted line) with the mass spectrum of 
protogalactic fragments inferred from the 28 globular cluster systems studied here (filled circles). Following Moore et al. (1999), we express 
the mass function of dark matter halos as n{M) oc M.~^ . 



